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;_ " " Abstract 

^^ ' We discuss the SM Higgs discovery potential of the LHC in the channel pp — > i/ + jet ^ 77 + jet 

0^ ■ when the jet is observed at sufficiently high pt and a small rapidity to be reliably identified. We 

calculate all the signal subprocesses and the irreducible background with realistic kinematical cuts. 
^ \ The reducible QCD background is also estimated. 

li ' We conclude that the channel 77 + jet can give about 120-200 signal events for Higgs mass 

[] ' Mh =100-140 GeV at the integrated luminosity of 30fb^^. This signal rate should be compared 

with only '~330-600 events for the irreducible background per two-photon invariant mass interval of 
2 GeV. We estimate the QCD reducible background at the level of < 20% of the irreducible one. 
Thus, one may hope that the Higgs boson can be discovered already during the LHC operation at a 
low luminosity. At a high luminosity of lO^'' cm~^s~^ the observation of several hundreds of high pt 
^ ' Higgs bosons in this channel will be possible with significance higher than 15 for L ~ lOOfb"^. 
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It is well known that the observation of a hght Higgs boson [Mh < 140 GeV) at the LHC collider in the 

inclusive channel pp ^- 77 -I- X is not easy 0, 0, 0] . Even if we can reduce the misidentification of two 

jets as photons to a small level, the irreducible background from the qq -^ 77 and gg -^ 77 subprocesses 

1-5 I rapidly increases for smaller 77 pair invariant masses, and it is necessary to separate a rather elusive 

Higgs boson signal from it. In this situation it is important to understand whether we can observe at the 

Q^l LHC any other Higgs boson production mechanisms. In this paper we are considering the Higgs discovery 

4j ■ potential in the channel pp —^ 77-I- jei when the jet is observed at a sufficiently high transverse momentum 

and small rapidity to be reliably identified. In this channel the signal cross section is much smaller in 

comparison with the inclusive pp — > 77-t-X case, but at the same time the situation with the background 

is undoubtedly much better. Transition to larger pt allows us to reduce the overwhelming backgrounds 

of the completely inclusive channel mentioned above. We consider the scalar boson production in the 

Cd ' framework of Standard Model (SM) with one Higgs doublet. 

The idea to look for the Higgs boson signal associated with a high pt jet in the final state was considered 
in Pi, where the matrix elements of signal subprocesses gg^g + H,gq^q + H and qq ^ g + H were 
calculated analytically in the leading order a^. In |^ the corresponding subprocesses were calculated 
for the case of CP-odd Higgs boson production within the minimal supersymmetric extension of the 
Standard Model (MSSM). Recently (see Ref. Q) the matrix elements for subprocesses gq -^ H^q' , H^q 
were calculated analytically for charged and neutral scalars in the framework of MSSM. However, in M, M 
decay channels of the Higgs were not considered and in [p| only the final state t^t~ + jet was discussed. 
In M the SM subprocesses were considered for the case of heavy Higgs boson decaying into the WW 
or ZZ pairs. Very promising numbers have been obtained recently in M for 77 -I- 2jets final state with 
two very forward jets (|?7jet| < 5), when Higgs scalars are produced by weak boson fusion mechanism. 
In this reaction 10-20 Higgs signal events could be observed with the significance equal to 3.5-7 at the 
integrated luminosity 10 fb^^. In |1^ (see also 0) the final state 77 -I- (> 2 jets) was simulated by means 
of PYTHIA generator [|llj in the realistic CMS detector environment. The events with relatively small 

*To appear in: Proceedings of Xllth International Workshop: High Energy Physics and Quantum Field Theory 
(QFTHEP-97), Samara, Russia, September 4-10, 1997. Preliminary version of this paper was included in [hi, where 
we studied also the possibilities of Higgs signal observation in the final state 77 -|- lepton for the LHC high luminosity 
regime. 



rapidity jets and the events with very forward jets, produced by weak boson fusion mechanism, were 
studied there (as we mentioned above, recent independent analysis of the fusion can be found in 0). It 
was shown that this signature has good prospects for the hght Higgs boson search giving ^^100-150 signal 
events at the integrated luminosity of 160 fb~^ (about one and a half year of LHC operation at the high 
luminosity of 10^'*cm~^s~"'^) with approximately the same number of background events. 

The final state 77 + jet, with high pt and small rapidity jet recoiling against the Higgs boson, has 
not been analysed and we discuss it in detail in Section 4. The possibility of charged particle tracks 
reconstruction in the central detector is an advantage of this signature {Ef'^ > 40 GeV and \rijet\ < 2.4). 
Final state jet reconstruction allows to determine more precisely the position of interaction vertex, thus 
giving the possibilities to improve two photon invariant mass resolution. This point could be important 
for the separation of background. We calculate all the signal subprocesses and the irreducible background 
applying realistic kinematical cuts. The reducible QCD background is estimated as well. 

We calculated cross sections and distributions with the help of CompHEP package jT2| |l^ . Exact one- 
loop matrix elements for the signal QCD subprocesses were implemented in the CompHEP FORTRAN 
output, thus we got a possibility to use CompHEP numerical module for integration over the phase 
space. Then, the vertex ggH in the signal QCD diagrams, including quark loops, was implemented in 
CompHEP as an effective point-like vertex together with the gggH vertex to ensure gauge inivariance of 
the amplitude. Thus we tested an accuracy of the effective Lagrangian approximation. 

The CompHEP package includes the code of adaptive MC integrator VEGAS |jlj] . For parton distri- 
bution functions we used the paramctrizations CTEQSm ||l||, CTEQ4m and CTEQ41 |l|], MRS-A' ||l|| 
and GRV92 (HO) |18]. The corresponding codes are also implemented in CompHEP. 

The paper is organized as follows. In Section 2 we list the values of physical parameters and some 
physical conventions. 

In Section 3 we discuss basic formulas for the H -^ 77 branching and the ggH vertex. To test the 
code written by using these basic formulas we compare results obtained by means of CompHEP, PYTHIA 
6.1 |ll[ and the programs HDECAY/HIGLU |jl9] for 7? ^- 77 and for the reaction pp ^ H ^ 77. 

In Section 4 we analyse the reaction pp — > 77 -I- jet. In Section 4.1 the signal QCD subprocesses 
gg —> H + g, gq —f H + q and qq —^ H + g are discussed. The corresponding formulas for their 
matrix elements in the limit of heavy i-quark mass are given, what turns out to be a precise enough 
approximation. Detailed analysis of the main signal and background subprocesses, g.g ^ iJ + 5 — > 77 + <? 
and gq ^ j + j + q, is presented in Section 4.2. We introduce different cuts looking for the possibilities to 
get a better significance of the Higgs signal in the 77-f jei final state. In Section 4.3 the contributions of 
all signal and irreducible background subprocesses are discussed. Then, in Section 4.4, we estimate the 
QCD reducible background. In Section 4.5 a general discussion for the reaction pp -^ 77 -I- jet is given. 

In Conclusion the prospects to search for the light Higgs boson at low LHC luminosity are summarised. 

2 Input parameters and physical conventions 

In our calculations we used y/s — 14 TeV and Mz = 91.1884 GeV, sin^ 0^, — 0.2237. For heavy quark 
masses we took the values nis = 0.2 GeV, TOc = 1.3 GeV, nib = 4.3 GeV and rrit = 175 GeV Effl R 

Most of the hard subprocesses that we are discussing involve QCD vertices. For the strong coupling 
constant we are using the normalization 

a,(Mz)= 0.118, 
where as runs in correspondence with the 2nd order QCD (NLO) formulas with flavour matching (see, 
e.g., [|0|, p.77). This normahzation assumes that A^^^^ = 226 MeV. For instance, as(120GeV) = 0.1133. 

In electroweak vertices we used a{Mz) — 1/128.9 [Q calculating electroweak couplings and the 
Fermi constant by means of the corresponding formulas with a{Mz), Mz and sin 9^ as independent 
parameters. However, it is more accurate to use a{me) = 1/137.036 in each vertex that includes an 
on-shell photon because of the absence of photon wave function renormalization in this case. Thus we 
introduce the correction factor (128.9/137.036)^ for each process with two photons in the final state. 

^We shall not discuss the sensitivity of our results to the value of t-quark mass. Recent CDF/DO data gives mt = 
175.6 it 5.5 GeV pi|. Rather small uncertainty of mt does not change significantly cross sections of the processes under 
consideration. 



3 Higgs decays and the effective Lagrangian 

It is well known that in the case when Mh < 140 GeV the Higgs boson decays dominantly into a bb pair. 
However, the QCD corrections to this decay partial width are large. So, for reactions where the ff ^ 77 
decay occurs it is important to take into account carefully these corrections because they increase the 
two-photon branching approximately by a factor of 2. In this section we also discuss all other decay 
channels related to the i/ — > 77 branching and the effective 77_ff Lagrangian which can be used for 
calculation of the partial width T^^^^^ . Similar effective ggH Lagrangian can be used in calculation of 
the partial width of Higgs decay into a gluon pair and the cross sections of signal QCD processes. 

In order to take into account the QCD corrections the NNLO formulas improved by the renormal- 
ization group method [ p2| were proposed for the running quark masses. We implemented in CompHEP 
the corresponding NLO formulas for the Higgs decay into quarks. This approximation gives results for 
partial widths higher than the NNLO formulas by ^ 10% and the corresponding correction factor has 
been introduced in our computer code. 

Electroweak corrections to the Higgs decays into quarks and leptons are small [g3| and we neglect 
them. 

Analytical formulas for the partial widths Tj^^^^-y and T^^gg in the leading one-loop approximation 
are known for many years p4| . They can be obtained with the help of the following effective Lagrangians: 

^1"h = K^Gl^G'^^'^H . (1) 

Here F^i, and Gf^jj are photon and gluon field strengths, correspondingly, and K is a factor accounting 
for high order QCD corrections. Analytical formulas for the effective coupling constants X^^h and \ggH 
in the leading order approximation can be found in Refs. [^^. The corresponding partial widths are 

\2 \2 

i H^77 - -^ ■ Mh , i H^gg - J^ ' -^^ ' M^ ■ 

Total QCD radiative corrections to the H ^ gg partial width can be collected in a factor with value 
'^ 1.7 |2^, |26, 27, ^. Note that the QCD radiative corrections to the partial width of i? ^ 77 are very 
small (less than 1%, see, e.g. |^ and references therein). 

Electroweak corrections to the partial widths of iJ — s- 77 and H ^ gg are small and can be neglected 

In the limit rrit -^ 00 and without contributions from the light quark loops the ggH effective coupling 
constant has a very simple form 

mt^oo 127rAf^y 

where the electroweak constant is a^ = g^/^^ = Oi{Mz)/ sin^ 9w As we shall see the effective Lagrangian 
(|l|) with the coupling constant (g) can be used in calculations of the signal subprocesses with acceptable 
accuracy. 

We remark that for a precise calculation of the _ff ^ 77 branching in the Higgs mass range considered 
one has to account also for the partial width Vh^ww , where W* means the off-shell M^-boson. At the 
mass values close to Mh — 140 GeV this branching and the branching to a bb pair are of the same order. 
The corresponding analytical formulas can be found in | pl| . 

All formulas for the partial widths mentioned above and the higher order QCD correction factors were 
implemented in the CompHEP codes. We checked the numbers obtained for the branching 7J ^ 77 by 
means of CompHEP, PYTHIA 6.1 and the program HDECAY ^■. 

Br{H -^ 77) • 10-3 



Mh GeV CompHEP HDECAY PYTHIA 



100 
120 

140 



1.999 


1.961 


2.065 


2.723 


2.707 


2.857 


2.148 


2.151 


2.327 



One can see a good agreement of our numbers with the results of the HDECAY program (where all 
known formulas for higher order corrections are used). At the same time one can see that the PYTHIA 
results are 5-8% higher. 

To test our code written for the effective coupling constant XggH (in leading order) we compared also 
the numbers obtained for the pp -^ H ^ jj cross section by means of CompHEP, PYTHIA 6.1 and 
HIGLU/HDECAY programs: 

a*°*(5g^H^77),fb 



package 


PDF set 


Mh = 100 GeV 


120 GeV 


140 GeV 


CompHEP 


CTEQSm 
GRV92(HO) 


37.18 
38.08 


41.03 
41.11 


26.85 
26.45 



HIGLU/HDECAY GRV92(HO) 39.59 40.78 25.17 

PYTHIA 6.1 CTEQ3m 38.96 40.97 26.65 

Here we used Q^ = M^ for the parton factorization scale. The strong coupling was taken fixed at 
the value of Higgs mass. The hard subprocess gg ^ H was calculated in leading order a^. Values of the 
H — > 77 branching used in these programs are given above. One can see that all three programs are in 
agreement with each other within 5%. 

4 Higgs boson signal in the reaction pp — > 77 + jet 

In this section we study systematically processes with the final state 77 + jet. 

First of all let us note that the signal subprocesses contributing to this reaction are QCD subprocesses, 
and it means that the cross sections should depend strongly on QCD parameters. There are three main 
sources of this dependence. The first one is defined by the evolution of the parton densities. Then, the 
hard subprocess cross sections depend on the running as- Finally, the iJ — > 77 branching depends on 
the Higgs total width, where the QCD corrections are large (see a short discussion in Section 3 and the 
corresponding references). In leading order all these three sources can be factorized. Moreover, due to 
a very small value of F^* in the Higgs mass range of 100-140 GeV the decay width of 7? — > 77 can be 
factorized, and the fixed value of strong coupling as{MH) can be used for evaluation of the H — > 77 
branching. However, it is well-known that for the reaction pp —^ H -^ '^j the dependence on the QCD 
normalization scale /i (used in calculations of the QCD corrections to hard subprocesses) and on the 
parton factorization scale Q is strong enough and next-to-leading analysis is necessary (see, e.g. [ p7| and 
references therein). For the reaction pp ^ H ^ 77 after including the NLO corrections, this theoretical 
uncertainty decreases considerably, showing only a '^ 15% remaining (fi, Q) sensitivity. Surely one can 
expect a similar effect also for the Higgs production at high pt- The self-consistent analysis requires 
the NLO corrections to hard subprocesses which are not known yet. Thus, today we cannot analyse 
the reaction pp — > 77 -(- jet in the complete NLO approximation. Therefore, in this paper some kind of 
combined accounting for the QCD effects is used: NLO PDF evolution and LO approximation for hard 
subprocesses (but with NLO running a^). 



4.1 Signal QCD subprocesses 

There are three hard QCD subprocesses with a signal from the Higgs boson contributing to the reaction 
PP ^' 11 -\- jet: 

gg^ H + g , gq ^ H + q , qq^ H + g. (3) 

In Fig. n^ Feynman diagrams contributing in leading order a^ are shown. Here the vertices ggH and 
gggH are substituted instead of the corresponding quark loops in exact calculations and represent the 
effective point-like vertices of the Lagrangian (|^J|) in the approximate calculations. 

The corresponding matrix elements were calculated analytically in |q|. In the limit nit — > 00 the 
differential cross sections are 
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dt 
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Here s, i and u are the standard Mandelstam variables for partons in their cm. system. These formulas 
can be derived from the effective Lagrangian (|l|,||). 

We found that the approximation (^ works well enough. For the channel gg this approximation gives 
numbers slightly smaller than the exact matrix elements - the difference is less than 7% at the transverse 
momentum cut pt > 40 GeV and \rijet\ < 2.4. For the channel gq this approximation overestimates the 
exact result by ~ 13% with the same cuts. For the channel qq the approximation (^ does not work at 
all, however, the contribution of this channel to the process pp —^ jj + jet is much smaller than the 
contribution of gg and gq channels and can be neglected. Thus, the precision of the effective Lagrangian 
approximation in calculation of the QCD signal processes (^ is better than 7%. This fact allows in 
principle to create a fast generator of signal events for the reaction pp — > 77 + jet. Detailed analysis of 
the approximation ^ will be published elsewhere p^ . 

We checked that the PYTHIA result is also in good agreement with exact calculations, giving cross 
section higher by ~ 6%. 

4.2 pp — > H -\- g — >■ 77 + g versus pp — >■ 7 + 7 + g 

In this section we consider the contribution of the signal subprocess gg —f H + g ^ 77 + g (four diagrams 
in Fig. la). We compare the corresponding cross section and distributions with those of the background 
subprocess gq —^ J + "f + q, q = u,d. Tree level Feynman diagrams contributing to this background are 
shown in Fig. ^. Let us stress that the two subprocesses considered give the main contributions to the 
signal and background. In this section we look for the possibilities to get a better significance of the 
Higgs signal varying the cuts on different variables. 

In the following we consider realistic cuts H, ^ for the final state 77 + jet (which we shall refer as to 
the (C) set of cuts) i^ 

(CI) two photons are required with pj > 40 GeV and [jy^j < 2.5 for each photon; 

(C2) photons are isolated from each other by Ai?(7i, 72) > 0.3; 

(C3) jet has high transverse energy E^"^ > 40 GeV and is centrally produced, \rijet\ < 2.4; 

(C4) jet is isolated from the photons by AR{jet,ji) > 0.3 and AR{jet,j2) > 0.3. 

In the following we shall vary the value of the kinematical cut on one of the variables starting from 
this basic set and leave a more tuned multivariable analysis for the future. 

It is clear that the signal significance depends strongly on the two-photon invariant mass resolution. 
This key instrumental characteristic at a low luminosity is expected to be ~ 1.3 GeV for the CMS electro- 
magnetic calorimeter [^ and ~ 3.1 GeV for the ATLAS calorimeter |j] {Mh = 100 GeV). These numbers 
are for the case of ~ 80% signal events reconstruction. Hereafter we shall integrate the background over 
77 invariant mass within the range Mh — AM^-y < Mj^ < Mh + AM^^ with AM^^ — 1 GcV. 

Cross sections of the subprocesses under discussion convoluted with various parametrizations of the 
parton distributions for the basic set of kinematical cuts (C) and for Mh = 120 GeV are 

CTEQ4m CTEQ3ni MRS-A' 



(S) , fb 4.32 4.65 4.44 

(B)/2GeV , fb 11.89 11.94 11.94 



Thereinafter we denote by (S) the signal process pp —> H + g —> 77 -I- g and by (B) the background 
process pp — > 7 -I- 7 -I- g. Here we use Q^ = Mjj + 2[El'^ Y ^-s the parton factorization scale. 

To estimate the NLO corrections we calculated these cross sections also with the leading order 
parametrization CTEQ41 and with strong coupling running in leading order. The LO results appear 
to be ~ 15% lower for the signal and only ~ 3% lower for the background. For the signal about 2/3 of 



^A_R = yj /^rp + A(/)2 is a separation between two particles in the azimuth angle - rapidity plane. 



this difference is due to the running of a^ , while for the background practically all the decrase is related 
to the PDF evolution. 

The parton cm. energy distributions presented in Fig. H show that almost all events have a small 
parton cm. energies vs < 300 GeV. So the phase space region for the process under discussion is 
characterized by the relations v^ ~ Mh ~ m-t- It means that we have events with very small values of 
the Bjorken variable x ^ i ■ \Q^^. For these small x values the recent DIS data from HERA |Q show 
some discrepancy with the existing parametrizations, especially for the gluon distribution. This data, as 
well as the recent DIS data from the NMC and E665 experiments and from precision measurements of 
the inclusive jet production at Tevatron, were used to improve the parton distributions by the CTEQ 
Collaboration [^. We have found that its latest parametrization CTEQ4m gives results lower than the 
previous set CTEQSm by ^ 7.5% for the signal and only by ~ 0.5% for the background subprocesses 
(see the Table above). We use the CTEQ4m parametrization {as{Mz) = 0.116, A^^^ = 202 MeV) in 
our calculations. Furthemore, we calculated the cross sections with different choices of the Q^ parameter 
which was used for both the parton factorization scale and the normalization for running Us (numbers in 
the Table below are for M// = 120 GeV): 



(S), fb 




4.32 


4.99 


5.73 


(B)/2GeV , 


fb 


11.89 


12.18 


12.68 



One can see a much stronger dependence on the choice of QCD parameter Q^ for the signal cross 
sections than for the background ones. Of course, it is caused by higher QCD order (a^) of the signal 
subprocesses, while the background is only of order of as- In the following we use Q^ = Mfj + 2{El'^ )^ 
which is intuitively more relevant to physics discussed here. 

Let us look at the rapidity interval between the jet and the Higgs boson due to which the minijet 
corrections p3| could give some enhancemenet factor if one considers the semi-inclusive process pp -^ 
H + jet + A^. In Fig. the distribution in the rapidity interval {yjet — Vh) is presented. One can see that 
this interval is less than ±3. So the corresponding minijet corrections should be small unlike in the case 
of SSC discussed in ||. 

Then, in Fig. one can see that varying the cuts on photon rapidities around set (C) values does 
not help to improve signal significance. However, a change of the jet rapidity acceptance does show some 
prospects (Fig. R). In Table 13(a) PI one can see that the significance increases when the jets with larger 
rapidities are detected. For example, we get a ^ 14% increase of significance if the \ri^\ < 2.4 cut is 
replaced by \ri^ < 4 with involving the very forward hadron calorimeter. 

The distribution in the jet transverse energy for the two processes discussed are shown in Fig. |8| for 
two values of the photon transverse momentum cut p] > 20 and 40 GeV. One can see that the application 
of a stronger cut improves the S/B ratio. This ratio increases also for a stronger Ej'' cut (see Fig. ^!q 
Let us see, however, what significance do we have. We have found that the variation of pi around set 
(C) values does not help to improve it. However, the application of a weaker cut on the analogous jet 
variable E^'' does improve the signal significance. For example, a replacement of the cut E'j^ > 40 GeV 
by Ef"^ > 30 GeV increases the significance by ~ 8% but leads to a small decrease of the ratio S/B, see 
Table |l](b). 

We analysed also a variety of other cuts, in particular, the cuts on transverse momenta of the photons, 
their relative angle in the two-photon rest frame and their separation cut. We treated also the cut on 
photon pair transverse momentum p]^ ^'^ , and different cuts related to the planarity features of events. 
Unfortunately these variables did not show any possibilities to get a better signal significance. 

4.3 Other signal and (irreducible) background processes 

A number of other subprocesses contribute to the reaction pp — > 77 + jet, both to the signal and 
background. We have calculated all of them within the framework described above. Results are collected 
in Table 0. 



^Here the effective significance ug/ ^/ob is presented rather than usually used Ng/y/Ng. To get the latter one should 
multiply the effective significance by the square root of the integrated luminosity. 

*Some suppression in the interval {Mu — E-j.'^ )/2 < pj < {Mjj + Bj*^ )/2 is connected with the cut applied on the jet 
transverse energy E^'^ > E^'^ . 



Signal subprocesses can be subdivided in two groups. In the first group there are QCD subprocesses (Bl) 
discussed in Section 4.1. NumericaUy the gq channels give about 12% of the main gg signal contribution, 
while the qq channels can be neglected. 

The second group of signal subprocesses is made up of electroweak reactions with Higgs production 
through WW or ZZ fusion, where one spectator quark produces the jet with high Ef. Of course, if one 
discusses the signature with only one jet one should assume that the set of kinematical cuts applied vetos 
one quark jet in these processes. Here one has also to account for a reaction where the Higgs boson is 
produced in association with the W- or Z-bosons decaying into quark pairs. In Fi g. p| Feynman diagrams 
for these subprocesses are shown. The W/Z fusion processes were calculated in Q, |35| , while the processes 
of W/Z associated production were considered in [g6| . In total the W/Z fusion and association rates are 
about 33-40% of the QCD signal channels for Mh = 100-140 GeV. Here we used Q^ ^ (My/2)2 as the 
parton factorization scale for the fusion processes {V denotes Vl^ or Z), and Q^ — {My + MuY for the 
associated production. 

Note that the process pp ^ H + t + i can also contribute to the final state 77 + jet. However, our 
preliminary estimate has shown that the corresponding rate is very small. So, to be conservative we do 
not take into account this signal channel in our analysis. 

In total, we obtain 4.1-6.6 fb total cross section of the signal subprocesses in the Higgs mass range 
Mh = 100-140 GeV. 

Let us now look at the irreducible background. As we already discussed in the previous section, the 
subprocess 59 ^ 7 + 7 + (? gives the main background. Another channel, qq ^ 7 + 7 + 5, gives about 
36% of the main gq contribution. Feynman diagrams of these subprocesses are shown in Fig. pi 

The one-loop background process gg ^ J + J + g exists, when the photons are radiated from the 
quark loop. Unfortunately this subprocess is not calculated yet. We estimate its cross section by means 
of PYTHIA 6.1 simulations. For this purpose one can initialize PYTHIA subprocess 114, gg ^ 7 + 7, 
and switch on gluon bremsstrahlung from the initial states. This is definitely only one of many other 
physical contributions to this process. The result of these simulations shows us that this background is 
about 2-4% of the overall contribution coming from the gq and qq channels. 

In total, the irreducible background amounts to 10.9-19.3 fb/2GeV for the Higgs mass range under 
discussion [Mh = 100 - 140 GeV). 

4.4 Reducible QCD background 

In this section wc discuss the results of PYTHIA 6.1 simulations for various QCD processes which could 
give a background due to radiation of photons from the fragmentating quarks or giuons. Note that the 
photon production mechanism is related in particular to a Tr^'-meson (and other neutrals) production - the 
energetic tt'^ 's could result in a photon pair which will be detected as a single photon in the electromagnetic 
calorimeter. Thus, an important characteristic is the 71*^/7 rejection factor of the calorimeter. For 
example, for CMS detector it is assumed P to be ~ 3 from the possibility to distinguish the two photon 
clusters. Nevertheless, we do not take into account this factor in our present analysis. 
One kind of reducible background is coming from the subprocesses 

gq-^i + g + q, 99 -^i + q + q, qq' ^ 1 + q{9) + (i'{9), 

in the cases when the final gluon or quark produces an energetic photon and the jet escapes the detection. 
The corresponding rates were estimated by the following algorithm. For definitness let us consider the 
subprocess gq ^ j + 9 + q- Firstly, we calculate its cross section by means of CompHEP imposing the 
kinematical cuts for a gluon the same as for a photon. We apply the set (C) of kinematical cuts on this 
3-body final state and integrate in the invariant mass of photon and gluon (regarded as a photon) over 
the interval of 2 GeV around the central point M^^ — Mh- The cross section turns out to be ct '^ 2 pb 
(cr' ~ 2.7 pb when the final quark is considered as a photon). Then the probability to get a high energy 
(> 40 GeV) photon from the fragmentating quark or gluon without further generation of a detectable jet 
is estimated as P{'^/qveto) ^ 2 • 10^^ and P{'^/gveto) ^ 0.3 • lO^^'Q. As a result we have an estimate of 
this background as [a x P(^/qveto) + cr' x P{'^/ gveto)] ~1.1 fb/2GeV. Corresponding contributions from 
the gg and qq' subprocesses are estimated in the same manner and the result is ^0.5 fb/2GcV for each 
channel. Altogether these channels give a background at the level of ^ 2.1 fb/2GeV. 



Another kind of reducible background could come from the subprocesses 

gq^-f + q, qq-^-f + g, 

when the second photon is produced during the quark or gluon fragmentation but this jet is still detected. 
Let us consider the gq channel. As a first step we evaluate with the help of CompHEP the cross section 
of this process, applying the pt > 40 GeV, \r]^\ < 2.5 and \r]q\ < 2.5 cuts, and it turns out to be of order 
of 2800 pb. For further PYTHIA simulations we switch on a fragmentation in this process and look for 
two photons (which can be, in particular, a photon plus tt*^) satisfying the corresponding cuts from the 
set (C). We require also that the invariant mass of this pair should be in the interval 90 < M^~^ < 140 
GeV. In addition we require these photons to be separated: no charged particles in the cone AR < 0.3 
are produced and the total energy of neutral particles in the conical area 0.04 < AR < 0.3 is less than 
2 GeV. Then we look at these events and try to find, with the help of the jet finder built in PYTHIA, 
a jet having proper Et and separated from the generated 'good' photons. We generated in total more 
than 3 • 10^ events and no good events were found. As a result we estimate this background to be <0.4 
fb/2GeV. The analogous estimate for the channel qq is much smaller (<0.01 fb/2GeV). 

From this latter simulations one can extract a parameter which is useful for estimate of the next 
type of reducible background. The probability to get an energetic photon from the fragmentating 40 
GeV jet, separated from each other so that the jet is still detectable, can be estimated at a level of 
Pil/j^tdet) ^ 4 • lO^** Q. So, this background from the subprocesses gq -^ "f + q can be estimated as 
2800pb X 4 • 10~^ X P(7, 7), where -P(7, 7) is the unknown probability for two photons (or the photon and 
tt") to be separated and to have their invariant mass within the interval Mn — lGeV < M^^ < Mn+lGeV. 
From the above simulations one obtains ^(7, 7) < 4 • 10^^. 

Third kind of reducible background could comes from the pure QCD subprocesses of 2 ^ 2 type, when 
both particles in the final state are gluons or quarks. We have to discuss these processes in connection 
with the probability to get two separated and energetic photons from the fragmentating quarks and 
gluons. There are possible contributions from the following subprocessesfl 

99-" 9{<i)+9{q), gq^g + q, qq' -* q{g) + q\g)- 

We estimate the corresponding rates by the following formula 



^ X P — ^ X P — ^ X P(7,7). 
\jetdetj Kjet^^toJ 

Here a stands for the cross section of the subprocess when the cuts pt > 40 GeV and |r^jet| < 2.5 are 
applied. For the gg channel a = 1.3 • 10^ pb. However, this enormous cross section is reduced to the 
insignificant level of < 1.4 • 10"'^ fb by probabilities to get the proper photons. For other channels a has 
values much smaller and we have rates of < 10~^ fb for the gq and of < 10~^ fb for the qq' ones. In total, 
this type of background gives a very small contribution of ~ 0.0025 fb/2GeV. 

Finally we consider reactions of the 2^3 type when all three particles in the final state are gluons or 
quarks. At the parton level these processes can be derived from the processes discussed in the previous 
paragraph taking into account gluon radiation. So, as in the previous case, there are three subgroups 
of such processes with the gg, gq and qq' initial states. Here we have to account for events where two 
separated and energetic photons are produced by fragmentating quarks or gluons which escape detection. 
Let us consider the subprocess gg —^ g + g + g- If we consider two of final gluons as photons and apply our 
basic set (C) of kinematical cuts the cross section is of order of 5.8 • 10'^ pb/2GeV. To get the necessary 
estimate we use the formula 

cr X P ( -^ — ] X P^ ^ 



yjet-uetoj XJ^tveto, 

In the gg case we have to multiply this formula by a combinatorial factor 3 due to three variants to 
produce two photons from different gluons. We neglect the possibility when two photons are produced 
from the same gluon jet while the other gluon jet escapes detection. So, this background is estimated 
at the level of ~ 0.016 fb/2GeV. The analogous estimates for other subprocesses of this type give 0.018 

^Note that there is some double counting for the first type of QCD reducible background analysed in the beginning 
of this section due to the photon radiation from the final quarks. We neglect this double counting in our rather rough 
estimates. The same remark is valid for the last type of the QCD background discussed in this section. 



fb/2GeV for gq and 0.002 fb/2GeV for the qq' channels. In total, this type of QCD background gives 
rather small contribution of ~ 0.036 fb/2GeV. 

Altogether four types of QCD reducible background processes contribute with a rate of < 2.2 fb/2GeV. 
This is not a negligible background and one has to take it into account analysing the reaction pp — > 

77 + jet. 

Now we can conclude that the QCD reducible background, connected with the probability to get 
a photon from the fragmentating gluons or quarks, is potentially dangerous, especially if we remember 
various uncertainties in our estimates. Nevertheless, this background turns out to be less than 20% of 
the irreducible background. 

4.5 Summary 

From Table one can see, first of all, that with an integrated luminosity of 10 fb^^ this channel can 
give 40-70 signal events with a number of background events only by a factor of 2-3 higher. This result 
was obtained for the basic set of kinematical cuts (C) (see Section 4.2) and for resolution determined 
mass interval Mh ± 1 GcV in the two-photon invariant mass M^^. The significances are Ns/^/Nb ~4.0, 
5.3 and 4.1 for Mh = 100, 120 and 140 GeV, respectively, showing good prospects for discovery of the 
Higgs boson in this reaction at low LHC luminosity, when total statistics collected is ~ 30 fb^^. These 
results also mean that each year of LHC operation at high luminosity of 10'^^cm~^s"^ will give hundreds 
of events with high pt Higgs bosons associated with high pt jet in central detector. The signal significance 
will be of order of 15 in this case. 

Most improtant advantage of the reaction pp —^ 77 -I- jet in comparison with the 'standard' channel 
PP ""*■ 77 is significant improvement of the Signal/Background ratio. Let us remind that in pp — > 77 
this ratio is ~ 1/15 (calculated in the LO approximation and with 2 GeV photon invariant mass intevral 
for the background, as in our case). For the reaction pp —f jj + jet the Signal/Background ratio is 
^ 1/2 — 1/3. Signal significance Ns/\/Nb is approximately the same for both channels. 

However, we should stress that in present analysis we did not account for various factors which can 
change the Higgs significance considerably. Let us review these factors and start from those which should 
help to improve the values obtained for the Higgs signal significance. 

First, the results presented in Table || for hard subprocesses have been obtained in leading orders 
al for the QCD signal processes and as for the background. As we have mentioned in Section 4.1 the 
next-to-leading QCD corrections are unknown for these subprocesses .R Notice that the NLO corrections 
to the hard subprocesses in the pp —^ jj + X inclusive reaction increase the Higgs production cross 
section by about 60% (see [M and the references therein). Of course, it is very probable that the NLO 
corrections enhance somehow the rates of background processes as well. However, due to a lower strong 
coupling order of the background processes one can assume their QCD corrections to be smaller than for 
the signal. 

There are possible ways to improve the significance of the 77 + jet channel by applying weaker jet 
cuts. In particular, if the detection of jet is possible at large rapidities, up to Irjjetl ~ 4 (involving the very 
forward hadron calorimeter), the significance will increase by ~ 14%. Furthermore, if the determination 
of jet is reliable with E^'^ > 30 GeV then the significance will increase by ~ 8% in comparison with the 
case of basic kinematical cuts (C). 

Let us remember now that for the CMS PbW04 electromagnetic calorimeter the two-photon invariant 
mass resolution is expected to be significantly better than 1 GeV in low luminosity running. For example, 
in p the following two-photon mass resolutions were considered for the reaction pp ^ 77 -f X at low 
luminosity: a^ = 540 and 600 MeV for Mh = 110 and 130 GeV, respectively. Surely, in this case the 
Higgs significance will increase noticeably compared to the estimates made earlier in this paper. 

However, there are also many factors reducing signal significance. One of them is that the one-loop 
subprocess gg — > 7 + 7 + 5 contributing to the background, unfortunately, is not still calculated - its 
exact matrix element is unknown. Our estimate of its rate is rather small, ~ 0.2-0.8 fb/2GeV. However, 
these numbers have been obtained by means of PYTHIA simulations for the process 35 ^ 7 + 7 where 
the final gluon is radiating from the initial state only. We expect the rate of this background to be higher 
due to other mechanisms of the final gluon radiation. 

^Some progress has appeared recently for the signal, see Rcf. Ist). 



Furthermore, the reducible QCD background is potentially dangerous as a rule at hadron colliders. 
We have estimated this background at the level of <2.2 fb/2GeV. Of course, more accurate evaluation of 
this background would be desirable taking into account the detailed simulation of the ATLAS and CMS 
detectors. Due to this discussion one can remember that similar analysis for the inclusive pp — > 77 + X 
reaction showed y, §| that the reducible QCD background can be suppressed at the level of < 10% of 
the irreducible bias. 

The effeciency of photons detection was not taken into account in our estimates. In our case the 
photons have large transverse momenta p] > 40 GeV and radiated with small rapidities \rij\ < 2.5. So 
one can hope that the photon detection efficiency will be high enough. At present time for both ATLAS 
and CMS detectors this efhciency is estimated at the level « 0.7 — 0.8. Experimental efficiency of hard 
jet detection was also not taken into account. 

Finally, combining all these factors one can hope that they could at least compensate each other, and 
our general conclusion about a possibility to discover the light Higgs boson during the LHC operation at 
low luminosity is still kept. 

Conclusions 

The channel 77 + jei gives very promising discovery possibilities for the Higgs boson with a mass of 100- 
140 GeV during the LHC operation at low luminosity of ^ lO'^'^cm^'^s^^. With an integrated luminosity 
of 30 fb~^ 120-200 signal events could be observed with 330-600 background events per 2 GeV interval 
of two photon invariant mass, showing a signal significance of '--^ 7 — 8.5. For high LHC luminosity of 
10'^'*cm~^s~^ hundreds of Higgs bosons will be observed with a significance of order of 15 for one year 
data accumulation. The reducible QCD background is at the level of <20% of the irreducible one. The 
next-to-leading QCD corrections, not calculated at present time, could enhance the signal more than 
the background. Higher two-photon invariant mass resolution that is now under discussion for the CMS 
PbW04 calorimeter (see ||], a^^ ^ 500 — 600 MeV at Mh = 110— 130 GeV) could increase the significance 
of signal. On the other hand, the efficiencies of photon and jet detection that we did not take into account 
decrease the signal significance. These factors (as well as others discussed above) work in the opposite 
directions and it is probable that they compensate each other. 
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Table 1: Cross sections of the main signal (S) and background (B) processes pp-^H + g^rjj + g 
and pp ^ J + "f + q, as a function of the jet a) rapidity cuts and b) transverse energy. The basic set 
(C) of cuts is imposed on other variables. Mh ~ 120 GeV. PDF set CTEQ4m is used. All subprocesses 
are evaluated in QCD leading order and the 2nd order running as with normalization as{Mz) = 0.118. 
Q^ = M|f + 2{El'^ Y is taken as the QCD scale for running as and as the parton factorization scale. The 
77 invariant mass for the background is integrated over the range Mh — AM^-y < M^^ < Mh + AM^^ 
with AM^^ = 1 GeV. 
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Hard subprocess 
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Table 2: Summary for the pp -^ j'j + jet reaction. Contributions of different subprocesses are shown 
for the basic set of kinematical cuts (C) (see Section 4.2). PDF set CTEQ4m is used. All subprocesses 
are evaluated in QCD leading order and the 2nd order running as with normalization as{Mz) = 0.118. 
Q^ = Mfj + 2{Ej'^*)'^ was taken as the QCD scale for running as and as the parton factorization scale 
in hard QCD subprocesses. The values Q^ = (My/2)2 and Q^ = (jVfy + Mh)^, where V = (W, Z), are 
used for the parton factorization scale in the W/Z fusion and W/Z associated Higgs signal procceses, 
correspondingly. The 77 invariant mass for the background is integrated over the range M^ — AM^^ < 

1 GeV. Asterisk '*' marks the results of the PYTHIA 6.1 simulations. 



M^^ < Mh + AM^^ with AM^^ 



Other results are obtained with the help of CompHEP package. 
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Figure 1: Feynman diagrams for the QCD signal subprocesses: a.) gg —^ H + g , h) gq —^ H + q and 
c) QQ ~^ H + g. In SM the ggH and gggH vertices include quark loops. In the effective Lagrangian 
approximation mQ) these vertices are point-like with the coupling constant \ggH- 
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Figure 2: Feynman diagrams for the signal subprocesses with a) WW and ZW fusion mechanisms of the 
Higgs boson production, b) associated ly-boson and c) associated Z-boson production. 
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Figure 3: Feynman diagrams for the subprocesses a) gq ^ j + j + q and b) QQ — > 7 + 7 + g contributing 
to the background. 



15 



da / dV^ fb/GeV 



u.^ 


- 




- 


(^) 


0.18 


— 


:• \^) 




- 


j| (B) 


0.16 


— 




0.14 


- 




0.12 


- 


'- 


0.1 


- 


■■ 1 

1 
1 


0.08 


- 


1 


0.06 


: 




0.04 


- 


1 
1 ., 

~i 
1 


0.02 


1 


"""^^^^^-^^^i-^^ 




III 


""^^ ' ' " ■=--^T^=-'" 1 - 1 ^__ 


i\ 


iiiiiiiiiiiiiiiiiiiiiiliiiiFilin- ■|'~'^r 1 -■ 1- H - 



100 



150 



200 



250 



300 350 



400 



450 



500 



550 



600 



Vh GeV 



Figure 4: Distributions in the parton cm. energy ys for the main signal (S) and background (B) processes 
pp -^ H+g -^ 77+5 and pp -^ j+j+q. Energy ^/s = 14 TcV and the Higgs boson mass Mh — 120 GeV. 
CTEQ4m PDF set and Q^ = M^ + 2{El''*'Y are used. The basic set of kinematical cuts (C) is imposed. 
The 77 invariant mass for the background is integrated over the range M/f — AM^-y < M^-^ < Mh + AM^^ 
with AM-y^ = 1 GeV. 
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Figure 5: Distribution in rapidity interval yjet — Vh for the main signal (S) process pp ^ H + g ^ 
77 + g. Energy y/s ^ 14 TeV and the Higgs boson mass Mh = 120 GeV. CTEQ4m PDF set and 
Q^ ~ Mfj + 2{E^^ )^ are used. The basic set of kinematical cuts (C) is imposed. 
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Figure 6: Distributions in photon rapidity for the main signal (S) and background (B) processes pp -^ 
H + g — > 77 + (? and pp ^ -f + j + q. Energy y/s — 14 TeV and the Higgs boson mass Mh = 120 
GeV. CTEQ4m PDF set and Q^ = Mfj + 2(Sf *)2 are used. The basic set of kinematical cuts (C) 
is imposed on other variables. The 77 invariant mass for the background is integrated over the range 



Mh 



AM^^ < Mj^ < Mh 



AM-y^ with AM^^ = 1 GeV. 
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Figure 7: Distributions in jet rapidity for the main signal (S) and background (B) processes pp -^ 
H + g -^ 77 + 5 and pp ^ j + j + q. Energy y/s ^ 14 TeV and the Higgs boson mass Mh = 120 
GeV. CTEQ4m PDF set and Q^ = Mfj + 2(Sf *)2 are used. The basic set of kinematical cuts (C) 
was imposed on other variables. The 77 invariant mass for the background is integrated over the range 
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AM^^ < M^^ < Mh 



AM^^ with AM^^ = 1 GeV. 
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Figure 8: Distributions in jet transverse energy for the main signal (S) and background (B) processes 
pp^^H + g^jj + g and pp ^ j + j + q. Energy ^/s = 14 TeV and the Higgs boson mass Mh =■ 120 
GeV. CTEQ4m PDF set and Q^ ^ ]^p^ + 2(Sf *)2 are used. The basic set of kinematical cuts (C) 
is imposed on other variables. The 77 invariant mass for the background is integrated over the range 



AM^^ < M^^ < Mh 



AM^-y with AM^-y = 1 GeV. 
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Figure 9: Distributions in photon transverse momentum for the main signal (S) and background (B) 
processes pp —^ H + g ^ 77 + g and pp —^ j + j + q. Energy ^/s = 14 TeV and the Higgs boson mass 
Mh = 120 GeV. CTEQ4m PDF set and Q'^ ^ Mfj + 2(£;f *)2 are used. The basic set of kinematical 



cuts (C) is imposed on other variables. The 77 invariant mass for the background is integrated over the 
range Mh 



AM^-y < M-y^ < Mh 



AM^^ with AM-y^ = 1 GeV. 
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